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A THEORY  OF  SCATTERING  BY  SINUSOIDAL  METAL  SURFACES 


1.  INTRODUCTION 

The  problem  of  scattering  by  periodically-corrugated  metal  surfaces  is 
analyzed  using  an  approach  which  extends  the  physical  optics  approximation 
into  a rigorous  theory  and  is  amenable  to  efficient  computer  evaluation. 

The  theory  was  derived  in  the  course  of  a study  of  the  much  discussed 
polarization  problem  of  microwave  landing  systems.  One  aspect  of  this 
problem  concerns  scattering  by  large  periodic  metal  surfaces  near  run- 
ways--such  as  hangar  doors --and  the  question  as  to  whether  the  specular 
reflection  coefficient  of  such  surfaces  can  be  significantly  reduced  by 
appropriate  choice  of  polarization. 

Numerical  data  obtained  by  computer  evaluation  of  the  theory  show  that 
horizontal  polarization  generally  leads  to  substantially  less  specular  reflec- 
tion than  does  vertical  polarization,  particularly  in  the  interesting  range 
of  low  angles  of  incidence,  i.  e.  , near  grazing.  (It  is  implied  here  that  the 
grooves  of  the  corrugated  surface  run  in  the  vertical  direction.  ) Further- 
more, it  was  found  that  circular  polarization  is  highly  effective  in  suppress- 
ing higher  order  grating  lobes  of  the  same  polarization  (circular  with  the 
same  sense  of  rotation)  as  the  incident  wave. 

A detailed  presentation  of  the  new  theory  will  be  given  in  a forthcoming 
technical  report.  For  a comprehensive  literature  survey  of  previous  work 
in  the  field  of  periodic  surface  scattering,  we  refer  to  a report  by  Tong  and 
Senior  [1J. 

2.  SPACE  HARMONICS  REPRESENTATION  OF  THE  SCATTER  FIELD 

Consider  a periodic  metal  surface  with  a sinusoidal  height  profile 

zo  = h sin  ( 27T  T ) (1) 


as  shown  in  Fig.  1.  A plane  wave  of  dependence  exp(ia)t)  and  amplitude  E 
is  assumed  to  be  incident  from  the  direction  0 . For  the  TE -polarization 


[1J  T.  C-H.  Tong  and  T.  B.  A.  Senior,  "Scattering  of  electromagnetic 
waves  by  a periodic  surface  with  arbitrary  profile,  " The  University 
of  Michigan,  Scientific  Report  No.  13,  April  1972,  Prepared  for  Air 
Force  Cambridge  Research  Laboratories , Bedford,  Mass.  01730, 
under  Contract  No.  F19&28-68-C -0071. 
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Surface  grooves  parallel  to 


. Periodic  metal  surface  illuminated  by  plane  wave:  Coordinates 
and  geometrical  parameters  used  in  analysis  of  scatter  problem. 
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where  the  propagation  angles  0^  are  determined  by 


sine  = sine  + m -r 
m a 


(4a) 


Apart  from  the  common  phase  factor  exp  {-ikxsin  0),  all  of  these  plane 
waves  are  periodic  with  d,  the  groove  width  of  the  metal  surface  (Floquet's 
theorem);  they  have  therefore  been  called  "space  harmonics. 11  For 
m-values  in  the  range 

(4b) 


- — (1  + sine)  < m < + — (1  - sine), 

A A 


the  propagation  angles  0 are  real,  and  the  associated  space  harmonics 


are  of  the  propagating  type;  for  the  remaining  m-values,  0^  is  complex. 


and  the  space  harmonics  are  of  the  evanescent  type,  that  is,  they  decrease 
exponentially  in  the  positive  z -direction.  Clearly,  only  propagating  space 
harmonics  transmit  real  power,  while  evanescent  space  harmonics  exhibit 
purely  reactive  power.  Conservation  of  energy  requires  for  TE -polarization 
that 


m 


P P*  rT 

L .E  cos©  = T E *E*  cos0 
Jr  m ni  m 

1 


(5a) 


and  for  TM -polarization  that 


m. 


KP.HP*  cos 0 = 7 H * H*  c o s 0 
‘■mm  m 

ml 


(5b) 


Jr  / xr 

where  H = ve/p  E . The  summations  are  limited  to  the  propagating 
spectral  orders  defined  by  Eq.  (4b).  The  space  harmonic  of  spectral  order 
m - 0 propagates  in  the  direction  0^  = 0 and  hence  is  termed  the  "specular 


reflected  wave. 


The  complex  amplitudes  E and  H of  the  space  harmonics  can  be 


m 


expressed  in  terms  of  the  current  distribution  on  the  metal  surface; 


VF  Em  = ' H 1 W expl1k(x0  s1"em  * z0  cose„  )}  X (6a) 
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Hm  “ + 2d  I*0"*  7 Ix.z(xO)  exprik(x0  sinqm  + zQ  cosaj  1 x (6b) 


V 2 . 


x (1  - z'  tg  O(n)dx0 


where  z^  is  given  by  Eq.  (1)  and 

zo-^-2"  5 cos  <2”T>  * 


The  symbols  I and  I denote  the  current  densities  in  the  TE-  and  TM- 
y x,  z 

cases,  respectively.  The  integrations  extend  over  one  period  of  the  metal 
surface. 

If  this  surface  period  is  large  compared  to  a wavelength,  it  may  be 

assumed  that  the  current  densities  are  determined  primarily  by  local 

effects,  and  the  physical  optics  approximation  may  be  used.  In  other  words, 

I and  I at  each  surface  element  may  be  approximated  by  the  currents 
y x,  z 

which  would  flow  in  the  element  if  it  were  part  of  an  infinite  metal  plane 
tangent  to  the  surface  zQ  at  the  element  location.  Accordingly, 


*y  (V  ^ +2 


Hx^X0,20*  * 70  Hz^X0,Z0^ 


I + z 


p cosA  + z sinB 


1 + *' 


exp(-ik(xg  sinA  - z.  cosA)) 
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E exp(-ik(x 


sin9  - zQ  cosr)) 


(7b) 
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where  z^  and  z^  are  the  height  and  slope,  respectively,  of  the  surface 
profile  as  given  by  Eqs.  (1)  and  (6c). 

Inserting  (7a)  and  (7b)  into  (6a)  and  (6b),  respectively,  results  in 

explicit  expressions  for  the  scatter  field  amplitudes.  Shadowing  effects  can 

be  taken  into  account  by  using  the  current  densities  (7a)  and  (7b)  in  the 

directly  illuminated  surface  regions  only,  while  setting  I and  I to  zero 

y x,z 

in  the  shadow  regions.  Numerical  evaluations  show,  however,  that  the 
physical  optics  approximation  does  not  necessarily  lead  to  results  of  suf- 
ficient accuracy  even  when  the  surface  period  d is  of  the  order  of  several 
wavelengths.  Large  errors  occur,  particularly  when  the  direction  of 
arrival  approaches  grazing  incidence  ( 0 -+  900  ).  Hence,  to  obtain  reliable 
results,  rigorous  solution  of  the  scatter  problem  is  required. 


3.  SOLUTION  OF  SCATTER  PROBLEM  (OUTLINE  OF  NEW  APPROACH) 


Rigorous  solution  of  the  scatter  problem  requires  exact  calculation  of 

the  amplitudes  E and  H . Unfortunately,  the  space  harmonic  represen- 
m m 

tations  (3a)  and  (3b)  do  not  necessarily  converge  in  the  space  ranges  within 
the  grooves.  Hence,  in  general,  these  representations  cannot  be  used  to 
formulate  the  boundary  condition  at  the  metal  surface 


+ IS) 


tang 


0 


(8) 


which  would  provide  a straightforward  approach  to  solution  of  the  problem 
(Rayleigh's  method)  [2]  . 

In  recent  years,  several  rigorous  methods  have  been  developed  based 
on  alternate  representations  of  the  scatter  field,  i.  e.  , on  representations 
which  remain  valid  for  the  space  range  within  the  surface  grooves.  Typi- 
cally, these  methods  utilize  condition  (8)  and  lead  to  inhomogeneous 
Fredholm  integral  equations  for  the  current  distribution  on  the  metal 


[" 2 "]  For  a discussion  of  the  Rayleigh  method  and  its  range  of  validity, 
see  Tong  and  Senior  [1],  page  1. 
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surface  [3J  • A survey  of  these  methods  can  be  found  in  Ref.  [1]  . 

In  the  present  report,  a rigorous  method  is  described  which  does  not 

use  boundary  condition  (8)  explicitly.  Instead,  the  current  distribution  is 

found  by  formulating  a space  harmonic  representation  for  the  (hypothetical) 

field  radiated  into  the  half-space  below  the  current-carrying  surface  and 

then  postulating  that  this  field  cancels  the  primary  plane  wave  in  the  lower 

half-space.  This  approach,  which  is  related  to  a method  suggested  by 

Waterman  for  the  analysis  of  scattering  by  bodies  of  finite  dimensions  [ 4 J , 

is  outlined  in  more  detail  in  the  following.  The  first  step  is  to  expand  the 

current  densities  1 and  I into  the  series 
y x,  z 

y»o>  ■ K <*<>>  s •„  f„  <«o>  <9a) 

7 n 

Ix,z(x0>  " K.z  (x0}  £ bn  9n  • (9b) 

A A 

where  Iy  and  ^ are  the  physical  optics  approximations  (7a)  and  (7b), 

respectively,  or  suitable  modifications  of  these  approximations.  The 
functions  f^  and  g^  are  assumed  to  form  complete  systems  in  the  range 

d/2  — Xq  d/2.  To  satisfy  Floquet's  theorem  for  periodic  structures, 
they  must  be  periodic  in  d: 

'n  (*0  * d>  ’ W 

% <*0  * d>  * %(*0) 


for  all  n 


[5]  An  exception  is  an  elegant  recent  theory  by  Hessel  and  Shmoys,  wh^ch 
applies  to  periodic  surfaces  having  a rectangular  height  profile.  This 
theory  employs  a (waveguide)  mode  expansion  for  the  field  within  the 
surface  grooves.  A summary  of  this  theory  can  be  found  in  the 
following  report:  A.  Hessel  and  J.  Shmoys,  "Computer  analysis  of 
propagation/ reflection  phenomena,"  Polytechnic  Institute  of  Brooklyn 
(Farmingdale,  New  York),  Final  Report,  20  August  1973  under  Contract 
No.  DAAB07-73-M-2716,  Prepared  for  the  U.  S.  Army  Electronics 
Command,  Fort  Monmouth,  N.  J. 

[4]  P.  C.  Waterman,  "Matrix  formulation  of  electromagnetic  scattering," 
Proc.  IEEE,  vol.  53,  August  1965,  pp.  805-812. 
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a and  b are  (unknown)  coefficients, 
n n 

Suppose  these  coefficients  are  not  chosen  correctly.  Then  the  currents 
(9a)  and  (9b)  will  radiate  a field  not  only  into  the  half- space  above  the  per- 
iodic surface,  but  also  into  the  half-space  below.  (We  assume  temporarily 
that  the  metal  surface  is  removed,  but  that  the  current  distribution  is 
maintained  in  place  and  that  the  half-space  z < 0 has  the  same  free-space 
properties  as  the  half-space  z > 0).  Similar  to  the  field  in  the  upper  half- 
space, the  field  in  the  lower  half-space  can  be  represented  as  a superpo- 
sition of  space  harmonics  propagating  (or  evanescing)  away  from  the  sur- 
face. For  TE -polarization 


E exp (- i k (x  sin0  - zcos9  )) 

m mm" 


and  for  TM- polarization 

+ ’ 

H ■ L H exp(-ik(x  sinP  - zcos9  ))  . 
y m*-i-  m mm 

As  in  the  upper  half- space,  the  direction  angles  0 are  determined  by 

m ’ 

Eq,  (4a).  Hence,  the  propagation  directions  of  the  space  harmonics  in  the 

upper  and  lower  half -spaces  are  related  by  imaging  at  the  z-plane.  The 

complex  amplitudes  E and  H are  linear  functions  of  the  current 

m m 

coefficients  a^  and  b^,  respectively.  For  TE -polarization 


a U 
n mn 


- < m < + » 


U « - T77  — — ~\f^  f+2  I (x  ) f (x  ) ~\j 1 + z'  y 
mn  2d  cosO^  V^Jd  y0n0v  0 


v expf i k(xft  s i n0  - z cos9  ) )dxft  ; 

U m u mu 


(Ha) 


and  for  TM -polarization 


IL  ““  P 

T H - E y,  b V 
e m n mn 
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- «»  < m < + 


Oob) 


» i * m 0 ■ 


iW  J^2  (,0>  9n  ‘V  (l  * ‘0  '9  9J  x 

" 2 


x exof ik(x  sin0  - z.  cosfi  ) ldx„  . 
1 0 m 0 m 0 


(»lb) 


where  zQ  and  z^  are  given  by  Eqs.  (1)  and  (6c). 

The  correct  current  densities  I and  I , i.  e.  , those  flowing  in  a metal- 

y x,  z ° 

surface, are  now  found  by  postulating  that  the  field  in  the  lower  half-space 
is  identically  zero.  This  requires  that  the  zero-order  space  harmonic 
cancels  the  primary  wave  (both  waves  travel  in  the  same  direction)  and 
that  all  higher-order  space  harmonics  vanish.  Hence  for  TE -polarization 
r cp 

_ “E  for  m « 0 

Em  * j „ (12a) 

0 for  m -t  0 ; 

and  for  TM -polarization  [5  j 


•Vi 


for  m 


(12b) 


for  m ^ 0 . 


These  equations  in  combination  with  (10a)  and  (10b)  lead  to  linear  systems 

for  the  unknown  current  coefficients  a and  b : 

n n 


U a 
mn  n 


( 1 3a) 


72  v b 

n nrn  n 


(13b) 


for  - oo  5 m < + 


[ 5 J It  is  interesting  to  note  that  use  of  the  physical  optics  approximation, 
i.  e.  , the  approximate  current  densities  (7a)  and  (7b),  results  in  zero 

order  amplitudes  E^  = -E  and  = -^e/ y E^*,  which  satisfy  the 

zero  order  conditions  (12a)  and  (12b)  rigorously.  This  is  true  for  both 
versions  of  the  physical  optics  approximation,  i.  e.  , when  shadowing 
effects  are  taken  into  account,  and  when  they  are  neglected.  However, 
the  amplitudes  of  higher-order  space  harmonics  do  not  generally  vanish 
in  this  approximation.  9 
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where 
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1 for  m = 0 


L 0 for  m r 0 


It  is  expedient  to  choose  the  exponential  Fourier  system  for  the  functions 


f and  g 
n Bn 


fn(x0)  " 9n(x0)  = exp(’2  170  ~d) 

■ “ ^ n » integer  < + « 


(14) 


The  integrals  (11a)  and  (lib)  for  U and  V can  then  be  evaluated  in  closed 

mn  mn 

form 


C , ,.n+] 
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for  m ■ 0 


U * 

mn 
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(-D 


n+m+1  cosB 


J (6  ) for  m / 0 
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V = -i 
mn 


(15b) 
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1 - *7  (m  + n)  — J (3  ) for  m 4 0 

d cosO  - coso  nrt-n  “in' 


with  3 » 2tt  ^ (coso  - cosO  ) 

> m k m 


In  the  evaluation,  the  physical  optics  current  density  (7b)  was  used  for 


I and  a modified  version  of  (7a)  for  I : 
x,  z y 


Iy(x0)  " 2 


# 


coso 


V^; 


exp(-ik(XgSino  - z^coso))  . 


(16) 


10 


. J*  ’* 


^ tm>  i-#'*  . <l*  ■ ■“  <*#•*•  _t 


A 


^jp^rrswr 


It  has  been  found  that  this  modification  enhances  computational  accuracy. 


The  incorporation  of  expressions  (15a)  and  (15b)  for  U and  V into 

mn  mn 

the  linear  systems  (13a)  and  (13b),  respectively,  permits  efficient  computer 

evaluation  of  the  unknown  Fourier  coefficients  a and  b . A matrix  inver- 

n n 

sion  (or  equivalent  procedure)  coupled  with  a routine  for  calculating  Bessel 
functions  of  complex  argument  are  required,  but  no  time-consuming  numer- 
ical integrations  are  needed. 


After  computing  the  current  coefficients  a and  b , calculation  of  the 

n n 

scatter  field  amplitudes  E^  and  H for  the  upper  half-space  is  straight- 
forward. Upon  inserting  Eqs.  (9a),  (9b),  and  (14)  into  Eqs.  (6a)  and  (6b), 
integrals  are  obtained  which  again  can  be  evaluated  in  closed  form.  In  the 
case  of  TE -polarization 


/?m"2nx  (cosO+cosftJ  . (,7c) 

If  the  incident  plane  wave  is  circularly  rather  than  linearly  polarized, 
the  space  harmonics  of  the  scatter  field,  in  general,  will  be  polarized 
elliptically.  In  this  case,  it  is  convenient  to  split  the  space  harmonics  into 
their  circularly  polarized  component  waves.  The  amplitude  of  the  compon- 
ent wave  rotating  in  the  same  sense  as  the  incident  wave  is 

E“*2  "J  • CM 

while  the  amplitude  of  the  component  wave  rotating  in  the  opposite  sense  is 

E.«0"  2 (Em  "Y'T  Hm)  ’ <18b) 

where  E^  and  Hm  are  the  scatter  amplitudes  for  TE-  and  TM -polarization 
given  by  Eqs.  (17a)  and  (17b),  respectively. 
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4.  NUMERICAL  EVALUATION 


A computer  program  for  numerical  evaluation  of  the  theory  has  been 
written  L^j;  evaluations  have  been  performed  on  a Burroughs  B-5500 
Computer. 

Numerical  data  was  obtained  for  surface  profiles  of  actual  hangar  doors. 
Full-scale  traces  of  three  such  profiles  are  shown  in  Fig.  2.  It  is  apparent 
that  these  profiles  can  be  closely  approximated  by  sinusoids,  where  only 
Profile  3 requires  a moderate  correction,  as  indicated  by  the  dashed  lines. 
In  Table  1,  the  parameters  d and  h of  these  profiles  are  listed  normalized 
relative  to  a wavelength  of  X = 6 cm  (C-band). 


Table  1.  Period  d and  Amplitude  h 
of  Surface  Profiles  Shown  in  Fig.  2; 
d and  h are  normalized  relative  to 
a wavelength  * = 6 cm  (C-band). 


Profile  No. 

Surface  Parameters 

d/\ 

h/A 

1 

1.  3 

0. 1333 

2 

2.  5 

0.  375 

3 

2.  25 

0.417 

Figure  3 shows  the  scatter  angles  0 of  surface  Profile  1 as  a function 

m 

of  incidence  angle  0 [7]  . In  Figures  4a  through  7b,  the  powers  P^  of 

the  (propagating)  space  harmonics  of  this  surface  are  plotted  vs.  0 for 
TE-,  TM-,  and  circular  polarization  of  the  incident  plane  wave.  Figures 
8 through  12b  and  13  through  17b  show  the  corresponding  data  for  Profiles 
2 and  3.  P^  is  the  power  transmitted  by  the  m1-*1  space  harmonic  through 

unit  area  of  a plane  z = const,  supposing  that  the  incident  power  per  unit 
area  is  unity.  Thus 


[6]  The  computer  program  is  based  on  a slightly  modified  theory  using 

sine  and  cosine  functions,  rather  than  exponential  Fourier  functions 

for  f and  g . 
n  6 *  8n 

[ J Note  that  the  scatter  angles  0^  are  counted  positive  in  the  first 

quadrant  of  the  x,  z -plane,  while  the  incidence  angle  0 is  counted 
positive  in  the  second  quadrant  (see  Fig.  1). 
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SCATTER  POWER,  Pm  SCATTER  POWER 
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Fig.  4a.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  8 
of  incident  wave  (Profile  1;  Spectral  orders:  m * *1,  0). 
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INCIDENCE  ANGLE,  6 

Fig.  4b.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  8 
of  incident  wave  ( rofile  1;  Spectral  orders:  m - -1,  -2). 
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CIRCULAR  POLARIZATION 
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X 


INCIDENCE  ANGLE, 


Fig.  6a.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  6 

of  incident  wave  (Profile  1;  Spectral  orders:  m - *1,  0; 

Polarization:  Circular,  with  same  sense  of  rotation  as  incident  wave.), 


CIRCULAR  POLARIZAT I ON 


i - 1.3 
X 


2-0. 1333 

X 


0 | 30*  J 60° 

13.34*  32.58" 

INCIDENCE  ANGLE,  8 

Fig.  6b.  Powers  P^  of  propagating  space  harmonics  vs.  arrival  angle  6 

of  incident  wave  (Profile  1;  Spectral  orders:  m - -1,  -2; 

Polarization:  Circulai , with  same  sense  of  rotation  as  incident  wave). 
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INCIDENCE  ANGLE,  6 


Powers  Pm  of  propagating  space  harmonics  vs. 
of  incident  wave  (Profile  2;  Spectral  orders 


TE- POLARIZATION 


INCIDENCE  ANGLE 


Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  9 
of  incident  wave  (Profile  2;  Spectral  orders:  m = -1  through  -4) 
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30®  f 

36.87° 
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Fig.  11a.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  6 

of  incident  wave  (Profile  2;  Spectral  orders:  m = *2,  +1,  0; 
Polarization:  Circular,  with  same  sense  of  rotation  as  incid 


dent  wave) 


CIRCULAR  POLARIZATION 
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0 f 30®  | 60® 

11.54°  36.87® 

INCIDENCE  ANGLE,  9 

Fig.  lib.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  0 

of  incident  wave  (Profile  2.;  Spectral  orders:  m ■ -1  through  -4; 
Polarization:  Circular,  with  same  sense  of  rotation  as  incident  wave) 
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SCATTER  POWER 


CIRCULAR  POLARIZATION 


INCIDENCE  ANGLE,  6 


Fig.  12a.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  9 of  incident  wave 

(Profile  2;  Spectral  orders:  m = +2,  +1,  0;  Polarization:  Circular,  with  sense 
of  rotation  opposite  to  that  of  incident  wave). 


Fig.  12b.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  9 of  incident  wave 

(Profile  2;  Spectral  orders:  m ■ -1  through  -4;  Polarization:  Circular,  with 
sense  of  rotation  opposite  to  that  of  incident  wave). 
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TE- POLARIZATION 


INCIDENCE  ANGLE 


Powers  of  propagating  space  harmonics  vs.  arrival  angle  9 
of  incident  wave  (Profile  3;  Spectral  orders:  m * *2,  *1,  0, 


TE- POLARIZATION 


INCIDENCE  ANGLE 


Powers  Pm  of  propagating  space  harmonics  vs. 
of  incident  wave  (Profile  3;  Spectral  orders 
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TM-POLARIZATION 

7 ■ 

- • 0.417 

X 

J 

INCIDENCE  ANCLE,  6 

Fig.  15a.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  6 
of  incident  wave  (Profile  3;  Spectral  orders:  m » +2,  +1,  0, 


TM- POLARIZATION 


INCIDENCE  ANGLE,  6 


Fig.  15b.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  6 
of  incident  wave  (Profile  3;  Spectral  orders:  m • -2,  -3,  -4 
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INCIDFNCF  ANGI.E.  8 


Fig.  16b.  Powers  Pm  of  propagating  space  harmonics  vs.  arrival  angle  6 

of  incident  wave  (Profile  d;  Spectral  orders:  m » -2,  -3,  -4; 
Polarization:  Circular,  with  same  sense  of  rotation  as  incident  wave). 


SCATTER  POWER,  PC0  SCATTER  POWER , 


CIRCULAR  POLARIZATION 


33.75°  : 

INCIDENCE  ANGLE,  e 


Powers  P^  of  propagating  space  harmonics  vs.  a 

(Profile  3;  Spectral  orders:  m * -2,  -3,  -4;  Pi 
of  rotation  opposite  to  that  of  incident  wave) 


— - 


E *E*  cos  9 
mm  m 

EP* E^*  cos6 


for  TE-polar izat ion 


(19a) 


H “H  cose 
mm  m 

« .P  CP* 

— E *E  cose 

u 


for  TM-polarizat ion 


(19b) 


In  the  case  of  circular  polarization,  the  power  of  each  space  harmonic 
is  split  into  two  (orthogonal)  contributions: 


,CS 


_CS  _CS* 

E *E  cos6 
mm  m 

_P  rP* 

E «E  cose 


(19c) 


and 


_C0  cC0* 

rn  E • E cos6 
-CO  m m m 

m = „P  CP*  - 
c «E  cos9 


09d) 


which  correspond  to  circularly  polarized  waves  having  the  same  and  the 
opposite  sense  of  rotation  as  the  incident  wave. 

It  is  apparent  from  the  vs.  6 curves  that  the  specular  reflection 

coefficient  (P^)  is  significantly  smaller  for  TM-  than  for  TE -polarization, 

particularly  in  the  range  of  large  incidence  angles  0 >60°.  Hence,  the 
choice  of  horizontal  polarization  over  vertical  polarization  for  microwave 
landing  systems  will  significantly  reduce  the  probability  of  false  guidance 
caused  by  specular  reflection  from  large  periodic  metal  surfaces  near 
runways.  (It  is  assumed  here  that  the  grooves  of  the  surfaces  run  in  the 
vertical  direction.  ) One  of  the  objectives  of  the  present  study  was  to  con- 
firm this  experimentally-observed  effect  theoretically.  However,  when 

is  small,  one  of  the  higher  order  space  harmonics  may  become  dominant 
and  carry  a major  portion  of  the  incident  power.  Observe,  for  example, 
spectral  orders  m = -1  for  Profile  1 (Fig.  5b)  and  m = -3  for  Profile  3. 

In  the  latter  case,  the  dominant  space  harmonic  is  more  or  less  back- 
scattered  and  therefore  not  likely  to  cause  false  guidance  by  microwave 
landing  systems;  see  Fig.  13.  At  angles  very  close  to  grazing  incidence, 

9 90  , the  specular  reflection  coefficient,  of  course,  approaches  unity 

for  TM-  as  well  as  for  TE- polarization. 


<£9 


i or  ail  three  profiles  considered  here,  circular  polarization  is  very 
effective  in  suppressing  higher  order  space  harmonics  having  the  same 
polarization  (same  sense  of  rotation)  as  the  incident  wave.  Higher  order 
space  harmonics  of  opposite  rotational  sense  are  present  with  significant 
amplitudes;  however,  they  could  not  be  received  by  an  airborne  antenna 
polarized  to  receive  the  primary  microwave  beam. 

All  power  curves  show  Rayleigh-type  Wood  anomalies.  Near  those  angles 
of  incidence  0 where  particular  space  harmonics  are  at  grazing  angles  (and 
thus  change  their  character  from  propagating  to  evanescent  or  vice  versa), 
the  power  distribution  over  the  remaining  propagating  space  harmonics  is 
reordered  rapidly,  which  results  in  steep  peaks  or  valleys  in  the  power 
curves  [8]  . These  anomalies  are  pronounced  in  the  case  of  TM-  and 
weaker  for  TE -polarization,  an  effect  discussed  in  detail  in  the  literature 
on  optical  gratings.  The  condition  for  Wood  anomalies  to  occur  is 

6 m = *90  ; hence  with  Eq.  (4) 

0 = s i n ' (t  I -m  7)  . 

n 

These  incidence  angles  are  indicated  in  the  figures  by  small  arrows.  Note 
that  both  anomalies  of  Profile  2 are  double  anomalies  (two  space  harmonics 
at  grazing  angles  simultaneously);  therefore  a particularly  strong  effect 
should  be  expected. 


5.  NUMERICAL  ACCURACY  AND  LIMITATIONS  OF  PROGRAM 

Two  criteria  were  used  to  check  the  accuracy  of  numerical  results: 
a power  criterion  and  a reciprocity  criterion. 

Conservation  of  power  requires  that 

m 


P = I 


"1 


for  TE-  and  TM -polarization,  and  that 

T ('PCS^PC0' 

— \ m m / 

m| 


(20a) 


(20b) 


[8J  The^curves  for  Profile  3 were  drawn  using  data  points  at  P intervals 
of  5 . Therefore,  a certain  amount  of  conjecturing  was  required 
near  anomalies,  and  the  curves  may  not  be  accurate.  In  plotting  the 
curves  for  Profiles  1 and  2,  data  points  at  0.  5°-intervals  were  used 
near  anomalies.  These  curves  should  be  reliable. 
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for  circular  polarization;  the  summations  extend  over  all  propagating  spec- 
tral orders.  In  general,  it  has  been  found  that  criterion  (ZOa)  is  satisfied 
with  higher  accuracy  for  TE-  than  for  TM -polarization;  but  exceptions  are 
not  uncommon.  The  relative  error,  i.  e.  , the  deviation  of  the  left-hand  side 
of  Eq.  (ZOa)  from  unity,  is  typically  in  the  order  of  1%  or  less  for  Profiles  Z 
and  3,  although  in  a few  cases  larger  errors  occur  for  TM -polarization. 

F or  Profile  1,  the  accuracy  is  very  high,  the  error  being  in  the  order  of 

10 '4  %. 


The  reciprocity  theorem  applied  to  the  problem  of  scattering  by  periodic 
surfaces  leads  to  the  following  observations.  Consider  two  plane  waves 
incident  from  directions  0 and  0^  . If  the  first  of  these  waves  scatters  a 
space  harmonic  (of  order  m)  into  the  direction  of  the  second,  such  that 
a b 

= - 0 , then  the  second  wave  generates  a space  harmonic  of  the  same 

spectral  order  m propagating  in  the  direction  of  the  first  incident  wave; 
b 2l 

hence  0^  = - 0 [9]  and  the  complex  amplitudes  of  the  two  space  harmonics 

are  related  by: 


ca  a h b 
E cosR  * E cosR 

mm  fn 


(21a) 


and 


..a  .a 

H cosfl 
m m 


m 


cos9 


(21b) 


where  the  first  equation  is  used  in  the  TE  case  and  the  second  in  the  TM  case. 
We  have  assumed  that  the  incident  waves  have  equal  amplitude  and  phase. 


Criteria  (ZOa)  and  (Zl)  supplement  each  other.  The  former  criterion 
checks  on  the  accumulated  power  of  all  propagating  space  harmonics,  the 
latter  criterion  provides  a check  on  amplitude  and  phase  of  individual  space 
harmonics.  In  general,  it  has  been  observed  that  when  the  power  criterion 
is  well  satisfied  (within  1%  or  better),  the  same  is  true  for  the  reciprocmy 
criterion. 

The  accuracy  of  the  numerical  results  has  been  found  to  be  critically 
dependent  on  the  groove  depth  Zh.  When  this  depth  is  less  than  one  wave- 
length, the  accuracy  as  measured  by  the  power  and  reciprocity  criteria  is, 
in  general,  acceptable;  when  Zh  <<  \ , the  accuracy  is  excellent.  However, 
when  Zh  approaches  or  exceeds  A , the  accuracy  sooner  or  later  breaks 
down  and  unacceptable  errors  result.  It  has  not  yet  been  determined  whether 


[9]  With  regard  to  the  minus  signs  in  the  relations  0 = - 0 and 

b a m 

6 = - 0 , see  footnote  [7]  on  page  1Z. 
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this  limitation  is  due  to  poor  convergence  of  the  linear  systems  (13a)  and 
(13b)  at  large  h,or  if  it  is  caused  by  subroutines  used  in  the  program  (such 
as  the  one  used  for  computing  Bessel  functions  of  complex  arguments). 


With  respect  to  groove  width  d,  no  limitation  of  the  program  has  been 
encountered.  Accuracy  tends  to  increase,  rather  than  decrease,  with  d. 
This  is  to  be  expected  since  the  underlying  analytical  method  is  basically 
an  extension  of  the  physical  optics  approximation. 


As  an  additional  check,  the  above  program  was  tested  against  an  exist- 
ing program  developed  by  Zaki  and  Neureuther  [10],  [11 J .Their  program 

is  based  on  a theory  which  uses  a Green's  function  to  derive  an  inhomogene- 
ous integral  equation  for  the  current  distribution  on  the  metal  surface. 

From  the  current  distribution,  the  amplitudes  of  the  space  harmonics  are 
computed.  This  integral  equation  is  solved  numerically,  assuming- -as  in 
the  present  theory- -a  sinusoidal  height  profile.  Because  of  high  computa- 
tional cost,  comparison  of  the  two  programs  had  to  be  limited  to  a small 
set  of  parameter  values,  namely, 


o = 0 and  60 


d/X  = 1.5,  2.5,  and  3.5 
h/\  = 0.  I and  0.2  5. 


v 

f 

ft 

. » 

* 

I 

r* 

£ 

s 


Hence,  the  comparison  must  be  considered  preliminary.  However,  the 
powers  P of  the  propagating  space  harmonics  calculated  with  these  pro- 
grams were  found  to  be  in  excellent  agreement  for  each  parameter  set. 

In  Table  2,  the  relative  error  according  to  the  power  criterion  (20b) 

e , |,  -EPm|  . ,00yo  , 
m 

is  given  together  with  the  computation  time  T for  each  run  [1-]  • For  the 
small  surface  depths  considered,  the  accuracy  of  the  new  program  (W-S), 


[10]  K.  A.  Zaki  and  A.  R.  Neureuther,  "Scattering  from  a perfectly 
conducting  surface  with  a sinusoidal  height  profile:  TE -polarization, " 
IEEE  Trans.  Antennas  & Progagation,  vol.  AP-19,  March  1971, 

pp.  208-214. 

[11]  , "Scattering  from  a perfectly  conducting  surface  with 
a sinusoidal  height  profile:  TM -polarization,"  IEEE  Trans.  Antennas 
&t  Propagation,  vol.  AP-19,  November  1971,  pp.  747-751. 

[12]  The  programs  were  run  on  a Burrough's  B-5500  Computer. 
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Table  2.  Comparison  of  Whitman-Schwering  (W-S)  Program  and  Zaki-Neureuther  (Z-N)  Program 
Concerning  Relative  Error  (according  to  power  criterion)  and  Computation  Time. 


Surface  Incidence 
Parameters  Angles 


d/A  | h/_X 


1.5  0.1 

2.5 

3.5 

1.5  0.1 

2.5 


1.5  0.25 

2.5 


(degrees) 

0 


Relative  Error  e in  Percent 


TE-Polari zation 


Computation  Time  T in  seconds 


0.73-10* 


0.44-10  0.15-10" 

0.40-10'4  0.33-10" 


TM-Polari zation 

TE-Polarization 

w-s 

Z-N 

W-S  Z-N 

0.77-10"® 

0.13 

193 

0.69-10"7 

10"2 

742 

0.27-10"8 

0.41-1 0*2 

1399 

0.46-10"4 

1.18 

188 

0.21-10"4 

0.78 

828 

0.29-10*® 

1.31 

1441 

0.30-10  0.26-10 


Subtotal  of  Computation  Time 


.77-10"2  0.12 

oe.in-2  1.42-10"® 


1.61-10"°  0.42-10" 


1.5  0.25 

2.5 


0.82-10"  1.72-10  0.55-10"‘ 

0.23  0.16-10'2  0.27 

3.07  0.34-10"2  4,14 


Subtotal  of  Computation  Time 


Total  Computation  Time 


"*.***■  rz 7t*r- 
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in  general,  exceeds  that  of  the  existing  program  (Z-N),  which,  however,  is 
fully  sufficient  for  all  practical  purposes.  With  increasing  h/  \ , the  rela- 

tive accuracy  of  the  two  programs  is  expected  to  invert:  the  Z-N  program  is 
probably  usable  at  h/  A values  substantially  in  excess  of  0.  5,  the  practical 
limit  of  the  W -S  program  at  least  in  its  present  form. 

The  computation  time  of  the  W -S  program  is  more  or  less  independent 
of  surface  parameters  and  incidence  angle;  therefore,  only  the  total  running 
time  is  presented.  The  average  time  per  run  is  £ 7 5 sec.  In  contrast,  the 
computation  time  of  the  Z-N  program  increases  substantially  with  d/A 
If  the  parameter  values  of  Table  Z are  used  as  a basis  for  comparison,  the 
W-S  program  runs,  on  the  average,  % 10  times  faster  than  the  Z-N  program. 
This  conclusion,  however,  should  be  regarded  with  caution,  since  it  is  based 
on  a small  sample  of  test  cases.  Moreover,  the  d/A  values  considered  here 
are  comparatively  large;  at  smaller  d/X  values  the  running  times  of  the  two 
programs  compare  more  favorably.  It  should  also  be  noted  that  the  Z-N 
program  includes  evaluation  of  the  physical  optics  approximation,  which  is 
not  considered  in  the  W-S  program. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 

(a)  A rigorous  theoretical  approach  to  the  problem  of  scattering  by 
periodic  metal  surfaces  has  been  discussed.  This  theory  reduces  the 
scatter  problem  to  the  familiar  task  of  solving  a system  of  linear  equations. 

In  the  case  considered  here,  the  periodic  surface  has  a sinusoidal  height 
profile  and  the  coefficients  of  the  linear  system  are  obtained  as  closed  form 
expressions  in  terms  of  Bessel  functions.  The  theory  is  amenable  to 
efficient  computer  evaluation. 

(b)  A second  class  of  periodic  surfaces  which  can  be  handled  efficiently 
by  the  new  approach  is  comprised  of  those  having  piecewise  linear  height 
profiles.  The  coefficient  matrices  of  the  corresponding  linear  systems 
reduce  to  closed  form  expressions  in  terms  of  exponential  functions.  A 
theory  for  trapezoidal  surfaces  is  already  available  and  has  been  programmed 
for  computer  evaluation  at  the  New  Jersey  Institute  of  Technology.  Evalua- 
tion is  in  progress. 

( c)  Numerical  evaluations  performed  for  sinusoidal  surfaces  with  per- 
iods of  the  order  of  a wavelength  show  that  the  specular  reflection  coefficient 
of  such  surfaces  is  significantly  smaller  for  TM-  than  for  TE -polarization, 

particularly  in  the  range  of  incidence  angles  0 > 60  . At  grazing  incidence, 

~t'  o 

0 = 90  , the  specular  reflection  coefficient,  of  course,  is  unity  for  both 
cases . 

(d)  Circular  polarization  was  found  to  be  highly  effective  in  suppressing 
higher  order  grating  lobes  of  the  same  polarization  as  the  incident  wave. 
Grating  lobes  circularly  polarized  with  a rotational  sense  opposite  to  that 
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of  the  incident  wave,  however,  are  present  with  significant  amplitudes. 


(e)  The  accuracy  of  the  W -S  computer  program  is  critically  dependent 
on  groove  depth.  Further  work  is  required  to  extend  the  applicability  of 
this  program  to  groove  depths  of  more  than  one  wavelength,  which  is  the 
present  limitation.  With  respect  to  groove  width,  no  limitation  of  the  pro- 
gram has  been  encountered  .(evaluations  have  been  performed  for  groove 
widths  up  to  10  A ). 

(f)  A systematic  comparison  of  the  new  program  with  the  Zaki  and 
Neureuther  program  and  other  existing  codes  (as  for  example,  the  Tong  and 
Senior  program)  is  recommended. 

(g)  Experiments  to  verify  theoretical  results  have  been  conducted  at 
ECOM.  These  experiments  were  designed  and  conducted  by  Dr.  J.  Mink 

[13]  . In  general,  good  agreement  between  theory  and  experiment  was 
obtained,  although  in  a few  cases --at  angles  of  incidence  near  grazing- - 
substantial  discrepancies  were  observed. 
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